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化 最 大 完工 时 间 为 目标 ， 设 计 出 一 种 具有 改进 解码 


方案 的 布谷 岛 算 法 对 其 进行 求解 。 ee 还 需要 兼顾 工人 的 加 工 情 况 ， 因 此 改 
进 了 传统 解码 方式 以 避免 机 器 和 工人 在 加 工时 间 上 的 冲突 ， 同 时 在 解码 时 尽 可 能 利用 机 器 和 工人 的 空闲 时 间 。 在 布 
We 将 布谷 鸟 种 群 随机 划分 为 3 个 子 群 ， 每 个 子 群 采用 不 同 Lévy 飞行 方式 独立 进行 寻 优 ， 并 通 
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分 算 子 实现 子 群 间 信息 交流 ， 不 仅 增强 了 算法 的 全 局 搜索 能 力也 平衡 了 算法 的 局 部 搜索 能 力 。 最 后 通过 基准 测 
A en 


验证 了 改进 布谷 鸟 算 法 和 改进 解码 方法 的 有 效 性 优越 性 。 
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Improved cuckoo algorithm for flexible job shop scheduling with dual resource constraints 
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China West Normal University, Nanchong Sichuan 637009, China) 


Abstract: Aiming at the flexible job shop scheduling problem with dual resource constraints (DRCFJSP) , this paper designs 


a cuckoo algorithm with an 1mproved decoding Scheme to solve it with the goal of optimizing the maximum completion time. 


Since DRCFJSP needs to consider the allocation of machines and the processing situation of workers, this paper improves the 
traditional decoding method to avoid the conflict of processing time between machines and workers, and makes use of the 
idle time of machines and workers as much as possible during decoding. Under the core framework of the cuckoo algorithm, 
this paper divides the cuckoo population into three subpopulations randomly, and each subpopulation adopts different Lévy 


flight methods to search for the optimum independently, and realizes the information exchange between subpopulations 
through the difference operator, which not only enhances the global search ability of the algorithm but also balances the local 
search ability of the algorithm. Through the experimental simulation analysis of the benchmark test case, the results show the 
effectiveness and superiority of the improved cuckoo algorithm and the improved decoding method. 
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双 资 源 约束 的 柔性 车 间 调 度 问题 (DRCFJSP) 


省 述 如 下 : 


ww 个 工人 W={W,WW,…W,} 操作 mw 台 机 器 M={M,M,,……,M,} 
对 个 工件 7=L72 7 进行 加 工 , 工件 6G=12… 站 包含 确 
定 的 有 道 工 序 {0,,0,.…,0.,}。 每 个 操作 的 加 工时 间 取 决 于 机 器 


和 工人 分 配 ， 


不 同 的 机 器 和 工人 可 能 会 导致 加 工时 间 不 同 。 
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数 为 
min Cu =max(CE,;) 
约束 条 件 如 下 : 
VieJ,VieJ],, > 2x =1 (1) 
keM weW 
VieJ,V ei,CE,, = SE + 2 2 (Ti) @Q) 
keM weW 
VieJ,veJ,Ch,;=SL,; + 2, > (Tye,) 0@3) 
keM weW 
VieJ,VvieJ{,2,...,n; —1),CE,; < SE jn (4) 
Vg,ie J,VheJ,,Viey,,vkeM, CE,, <SE,,+L(1— yw ) (5) 
Ve,ieJ,Vhe],, vie ,vweW,CE,, <SE,; +L(1- zw) (6) 
vkeM,S,>=0 (7) 
VkeM,vie], vie ,vweW,S: <SE +L(1— x,) (8) 
vkeM,VvieJ,vieJi,vweW,CE ;x < Fi (9) 
车 PW[O,)]=-1， 令 CEmmon=0， 则 有 VieJ.viey， 
SE =max(CEpyio, 1, CEpo, >, > 58 oo) (10) 
keM weW 
若 SW[O,j]=-1， 今 SLwio=Crw ， 则 有 VieJ,viey,, 
CL ; = max(SLsyto ,1 SLsto ,1 SLswto,]) (1 1) 
式 (1) 表 示 每 个 工序 只 能 在 一 台 机 器 上 由 一 个 工人 进行 


; 式 (2) 表 示 每 个 工序 的 最 早 完工 时 间 等 于 该 工序 最 早 开 
台 加 工时 间 与 所 需 加 工时 间 之 和 ; 式 (3) 表 示 最 晚 完 工时 间 等 
晚 开 始 加 工时 间 与 所 需 加 工时 间 之 和 ; 式 (4) 表 示 同 一 工 
件 上 工序 加 工 的 先后 顺序 约束 ; 式 (5) 表 示 每 个 机 器 任 一 时 刻 


最 多 


只 能 加 工 


道 工 序 ; 式 (6) 表 示 每 个 工人 任 一 时 刻 最 多 只 


能 加 工 一 道 工 序 ;， 式 (7) 表 示 所 有 机 器 不 能 提前 进行 工作 ， 从 


表 1 符号 变量 定义 0 时 刻 开 始 可 用 ; 式 (8) 工 件 在 机 器 上 进行 加 工时 机 器 必须 是 
Tab.1 Symbol variable definition 空闲 的 ; 式 (9) 机 器 必须 在 加 工 结束 后 才能 停止 ， 中 途 不 能 停 
符号 定义 止 ; 式 (10) 计 算 工序 最 早 开 始 加 工时 间 ; 式 (11) 计 算 工 序 最 晚 
n 工件 总 数 完工 时 间 。 
m 机 器 总 数 表 2 为 一 个 DRCFJSP 问题 实例 ， 该 实例 由 三 个 工件 
w 工人 总 数 个 机 器 和 两 个 工人 构成 。 表 中 给 出 了 每 道 工 序 在 某 个 机 器 上 
M 总 的 机 器 集 某 个 工人 加 工 所 需 的 时 间 ， 其 中 “-” 表 示 对 应 工人 无 法 操 
W 总 的 工人 集 作对 应 机 器 对 该 工序 进行 加 工 。 例 如 ，2. 无 法 由 工人 允 在 机 
0, 工件 i 的 第 j 道 工序 器 WM, 上 进行 加 工 ， 但 901 可 以 由 工人 W 在 M, 上 进行 加 工 且 
SE 工序 0, 最 早 开 始 加 工时 间 所 需 加 工时 间 为 1。 
CE,, 工序 0., 最 早 完工 时 间 表 2 DRCFJSP 问题 实例 
SL, 工序 0., 最 晚 开 始 加 工时 间 Tab.2 Examples of DRCFJSP problems 
Ch, 工序 O, 最 晚 完工 时 间 W, La 
S, 机 器 上 的 开始 加 工时 间 "02 Mm M M M M  M 
到 机 器 大 的 结束 加 工时 间 Oi 1 - 1 6 
六 一 个 足够 大 的 正 数 | 0 a s 2 六 
Ti 工序 0;j 在 机 器 上 由 工人 w 加 工 所 需 时 间 O, 5 - . E - 2 
pw 四 训 v 由 同一 工人 加 工 的 前 一 道 工 序 ， 著 v 为 该 工人 加 0 4 . 2 3 
的 第 一 道 工序 ， 则 PW[v]=-1 Os 4 6 - - 6 = 
5 加 训 v 由 同一 工人 加 工 的 后 一 道 工 序 ， 若 v 为 该 工人 加 A O,, 2 2 2 E 
的 最 后 一 道 工 序 ， 则 SW[v]=-1 Os 6 4 = = 4 1 
决策 变量 如 下 : ee 
不 如 果 工 序 Oy 在 机 器 k 上 由 工人 w 加 2 ”改进 布谷 鸟 算法 求解 DRCFJSP 
和 2.1 布谷 鸟 算法 
1 如 果 工 序 D,， 与 工序 0 都 在 机 器 k 上 加 工 ， 且 工序 布谷 鸟 算 法 作为 一 种 新 型 启发 式 搜 索 算 法 ， 主 要 基于 1 
| O,， 先 于 工序 0 个 更 新 策略 : 一 是 通过 Lévy 飞行 机 制 寻找 寄生 梨 穴 , 二 是 通 
0， 否 则 过 偏好 随机 游 走 策略 代替 被 发 现 的 鸟 梨 。 在 布谷 鸟 算法 中 每 
1， 如 果 工 序 D.， 与 工序 D.， 都 由 工人 中 加 工 ， 且 工 月 一 个 鸟 梨 代表 一 个 可 行 解 。 在 理想 状态 下 ， 布 谷 鸟 的 位 置 更 
Zehipw | 0。， 先 于 工序 0 新 公式 如 下 : 
0 否则 Xm = X!+0a®Levy(p) (12) 
以 最 小 化 最 大 完工 时 间 为 优化 目标 建立 数学 模型 ， 其 其 中 ，X! 表示 第 ii=12…,D) 个 鸟巢 在 第 : 代 的 位 置 ; @ 为 点 对 
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点 乘法 ; 4 通常 取 值 为 a=1 。 Zev(C) 为 随机 搜索 路 径 ， 服 从 
levy 分 布 : 
Lewy(B) ~ 1=17,1<B<3 (13) 
为 方便 运算 ， 文 献 [12] 使 用 式 (13) 产 生 /ev 随机 数 : 
Lew(p)=s*4 (14) 
lols 


其 中 ， KW 和。 均 服 从 标准 正 态 分 布 ，B=15。 


1 


| 0 
(下 于 的 xpx25) 
综合 上 述 公式 ， 布 谷 鸟 的 位 置 更 新 公式 如 下 : 
$xu 
Xr = XI! + x (X! — Xs) 
(16) 


按 发 现 概率 B 丢弃 部 分 解 后 ， 采 用 偏好 随机 游 走 重新 生 
成 相同 数量 的 新 解 ， 公 式 如 下 : 

X=X!+7y(X: — XE) (17) 
其 中 ，7Y 是 服从 均匀 分 布 的 缩放 因子 ，Xi,Xi 表示 第 1 代 的 两 

个 随机 数 。 
2.2 算法 编码 
于 标准 布谷 鸟 算法 适用 于 解决 连续 性 优化 问题 ， 无 法 
直接 将 其 用 于 解决 车 间 调 度 问题 ， 因 此 本 文 引入 最 小 位 置 规 
则 (smallest position value, SPV) 来 建立 个 体 与 实际 调度 的 映 
射 关 系 。 如 图 1 所 示 ，SPV 规则 是 将 原 向 量 进行 升序 排列 得 
到 新 向 量 ， 新 向 量 中 各 分 量 所 对 应 的 原 向 量 中 的 位 置 即 为 整 
数 编码 序列 。 本 文 在 算法 编码 时 先 随机 生成 工序 加 工序 列 
(OS), 然后 根据 所 生成 的 加 工序 列 ,对 机 器 和 工人 进行 分 配 。 
整个 算法 包含 了 三 层 序 列 ， 第 一 层 是 工序 加 工序 列 (0OS)， 第 
二 层 是 机 器 分 配 序列 (MA), 第 三 层 是 工人 分 配 序列 (WA)。 该 
编码 由 表 1 的 问题 实例 得 出 ， 表 示 3 个 工件 由 2 个 工人 在 3 
台 机 器 上 加 工 的 调度 ， 其 中 工件 1 包含 2 道 工序 、 工 件 2 包 
含 2 道 工序 ， 工 件 3 包含 3 道 工序 ， 工 件 的 加 工 顺序 为 : 
OO, 一 0 一 0 一 0 一 0 一 0 一 03 。 结 合 整个 序列 可 知 ， 工 序 
0 由 工人 2 在 机 器 3 上 加 工 ， 工 序 0;; 由 工人 2 在 机 器 3 上 加 
， 工 序 0 人 1 在 机 器 2 上 加 工 ， 其 余 工 序 依 此 类 推 。 
编码 带 | 1 I 2 2 3 3 3 


随机 值 | 0.21 | 0.37 -0.24|-0.55| 0.81 | 0.66 | 0.15 


| 4 本 ms 工人 
0S 2 2 3 1 1 号 3 
0, 1 CQ 2 Os 1 Ol OO CO， O33 
WA 2 这 1 1 1 2 2 
W WW W W W WW W 
MA 3 3 | 2 2 2 1 3 
Ms, M, M, M, M, M, M, 
图 1 算法 编码 过 程 
Fig.1 Algorithm coding process 


2.3 改进 解码 算法 

相 比 于 单 资源 约束 的 FJSP 问题 ，DRCFJSP 在 加 工时 不 
仅 受 前 一 道 工 序 的 结束 时 间 和 机 器 最 早 可 开始 加 工时 间 的 约 
束 ， 还 受到 工人 最 早 可 开始 加 工时 间 的 约束 。 因 此 本 文通 过 
对 插入 式 解 码 方案 进行 改进 ， 得 到 一 种 可 以 同时 对 机 器 和 工 
人 进行 主动 解码 的 改进 解码 方案 ， 改 进 解 码 方案 的 核心 是 有 
效 利 用 机 器 和 工人 加 工时 产生 的 空闲 时 间 ， 将 机 器 和 工人 提 
前 安排 到 合适 的 空闲 时 间 段 内 进行 加 工 ， 以 达到 缩短 整个 加 
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工时 间 的 目的 。 


记 工 件 i 的 第 j 道 工序 0; 的 开始 加 工时 间 为 
5 ， 结 束 加 工时 间 为 三 ， 该 工序 的 紧 前 完工 时 间 为 B17， 
工人 操作 机 器 上 对 工序 % 加工 所 需 时 间 为 7 。 在 解码 时 首 
先 需 要 判断 所 选 机 器 和 工人 的 加 工 情况 ， 再 根据 加 工 情况 判 
断 是 否 有 空闲 时 间 以 及 空闲 时 间 是 否 满足 插入 条 件 。 当 机 器 
和 工人 均 有 空闲 时 间 时 ， 有 以 下 三 种 可 能 出 现 的 情况 : 


情况 一 : 如 图 2 所 示 ， 机 器 和 工人 的 空闲 时 间 不 重合 ， 
对 此 无 法 进行 插入 操作 。 
工件 人 机 器 空闲 时 间 
工人 空闲 时 间 
a 
时 间 


图 2 第 一 种 情况 
Fig.2 The first situation 


情况 二 : 如 图 3 所 示 ， 机 器 和 工人 有 重合 的 空闲 时 间 段 
了 ,但 鳌 小 于 工序 所 需 加 工时 间 区， 即 五 <Zxe ， 因 此 无 法 进 
行 插入 操作 。 
工件 人 
机 器 空闲 时 间 
人 T, 工人 空闲 时 间 
SN ] 
时 间 
图 3 第 二 种 情况 
Fig. 3 The second Situation 
情况 三 : 如 图 4 所 示 ， 机 器 和 工人 有 重合 的 空闲 时 间 段 
， 且 五 大 于 等 于 工序 所 需 加 工时 间 区， 即 五 =me ， 因 此 可 
Ne。 
工件 人 
机 器 空闲 时 间 
T 工人 空闲 时 间 
时 间 


图 4 第 三 种 情况 
Fig.4 The third situation 

综合 解码 时 可 能 出 现 的 所 有 情况 ， 本 文 给 出 了 解码 算法 
的 伪 代 码 ， 如 算法 1 所 示 。 

算法 1 改进 解码 方法 伪 代 码 
输入 : 总 工序 数 TP ， 机 器 和 工人 可 加 工 信 息 。 
输出 : 工序 开始 时 间 SL 和 完工 时 间 E35, 机 器 加 工时 间 [SMW ,EX ] ， 
工人 加 工时 间 [SW ,EW.]。 
1:For w=1 to 7P 
2: 0;=0S(1p) ; fiag =0; % fiag 用 于 判断 是 否 插入 空 闪 时 间 
3: ”获取 0; 可 选 的 加 工 机 器 和 工人 s ， 加 工 所 需 时 间 Ti ; 
4: ”机 器 kk 和 工人 5 的 加 工时 间 [SMi,EMi]LSW,,EW,] 和 对 应 空闲 
时 间 [ASM', AEM' [ASW.,AEW.] ; 


5: If j=—1 % 工 件 第 一 道 工 序 

6: If [SA EM && [SW.,EW.]—= % 工 人 和 机 器 
都 还 未 进行 加 工 

7 S10; ED,=ST +tTank ;flag=1; 

8: End If 

9: Else 

16 : If [ASM,,AEM,]~= 儿 &K[ASW.,AEW.] 二 名 % 机 器 有 空 
闲 时 间 ， 工 人 无 空闲 时 间 

11: If ASM, >=max{EW,} && ASM': >= ET ;1 && ASMi + Ts <= AEM': 
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下 之 Ti=4SM ; ET; =S7+TDe ; flag =1 
13; End If 
14: Else If [4SM ,APENM 一 和 &&[ASW,AEW.]~= 名 % 机 


器 无 空闲 时 间 ， 工 人 有 空闲 时 间 


15 : If ASW, >=max{EM,} && ASW, >= ET ;, && ASW, +T;, <= AEW. 
16: STi=ASW, ; ET; = ST +Tns ; flag=1; 

17; End If 

18: Else If [ASM,,AEM.]~=® &&[ASW.,AEW]~=@ % 机 
器 和 工人 均 有 空闲 时 间 

19 : If max(ASM', ASW.) >= ET && max(ASM,, ASW.) + Ti <= 


min(AEM,, AEW.) 

20: ST; =max{ASM', ASW,} ; ET; =ST; +Ts ; flag =1; 
21: End If 

22 : Else % 机 器 和 工人 均 无 空闲 时 间 


23 : ST; =max{ET ;1, EMi, EW,} ; ET; =ST; +Tns ; flag =1; 
24: EndIf 

25: End1f 

26: If flag=0 

27: ST; =max{ET ;1, EMi, EW.} ; ET; =ST, + Ts ; flag =1; 
28: End1f 


29: ”更 新 工序 0; 的 最 早 完工 时 间 EZ 以 及 机 器 和 工人 的 加 工时 
间 和 空闲 时 间 ; 
30: End For 


2.4 算法 更 新 
2.4.1 改进 Lévy 飞行 


度 
索 
缺 
入 
进 
群 


计算 公式 如 下 


有 .OH 


较 差 的 两 个 个 体 。 
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改进 后 的 布谷 岛 算 法 如 算法 2 所 示 。 


算法 2 改进 布谷 鸟 算法 更 新 伪 代 码 


输入 : 迭代 次 数 iter ， 
当前 最 优 鸟 巢 信 息 Xevs 。 


输出 : 


: For i=]1 to iter 


to 3 


子 群 inest 。 


1 
2 For inest=1 
3 


End For 


每 个 子 群 采 
子 更 新 式 (12)(18)(19) 
4: 对 鸟巢 进行 评 
E 成 相同 数量 的 新 鸟 集 ; 


不 同步 长 因子 的 Lévy 飞行 进 和 


Xs 
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(20) 
表示 子 群 中 


了 更 新 ， 步 长 因 


价 并 按照 概率 R 舍弃 较 差 的 鸟巢 , 采 


区 人 人 


对 这 一 代 乌 站 
If 1/ 一 0 


和 
5 
6 
7 
8 
9 End If 


10: End For 
2.5 算法 流程 


Stepl 初始 化 参数 ， 布 谷 鸟 种 群 数 


iter ， 发 现 概率 PP 。 


进行 评价 ， 


更 新 当前 最 优 鸟巢 信 


按照 式 (28) 对 子 群 进行 差分 ; 


式 (17) 


息 和 


心 、 best,g 3» 


Step2 初始 化 种 群 ， 根 据 映射 规则 ， 将 岛 则 


n ， 最 大 迭代 次 数 


信息 转换 为 


包含 调度 信息 的 可 行 序列 ， 采 


在 标准 布谷 鸟 搜索 算法 中 , Lévy 飞行 的 步 长 越 长 搜索 精 
越 低 ， 适 用 于 全 局 探索 ， 步 长 越 短 搜索 精度 越 高 ， 局 部 搜 
能 力 越 强 。 由 于 标准 Lévy 飞行 的 步 长 因子 是 一 个 固定 值 ， 
乏 自 适应 性 ， 可 能 导致 算法 在 搜索 时 搜索 速度 慢 且 容易 陷 


子 群 。 


改进 解码 算法 对 


其 进行 解码 


得 到 完工 时 间 ， 保 留 当前 最 优 鸟巢 ， 


Step3 分 别 采 上 
行 公式 更 新 子 群 


局 部 最 优 。 针 对 这 种 情况 ， 本 文采 用 两 种 方法 对 步 长 因子 
行 改进 ， 通 过 将 布谷 鸟 种 群 随机 划分 为 三 个 子 群 ， 三 个 子 
分 别 采用 固定 步 长 因子 的 Lévy 飞行 方式 与 两 种 改进 自 适 


yy 


步 长 因子 的 Lévy 飞行 方式 进行 独立 的 更 新 操作 。 根 据 不 


同 
达 


~ 


优 


Lévy 飞行 在 当前 最 优 鸟巢 附近 进行 较为 


锁 


™ 


Xu 表示 当前 最 优 鸟巢 位 置 。 


进 


搜索 阶段 种 群 的 搜索 情况 ， 自 动 调整 步 长 因子 的 大 小 ， 以 
到 平衡 全 局 寻 优 速率 和 搜索 精度 的 关系 的 目的 。 

中 , 第 一 种 改进 方法 如 式 (18) 所 示 , 根据 当前 鸟巢 与 当前 最 
鸟巢 的 位 置 关 系 ， 对 步 长 因子 4 进行 自 适 应 调整 ， 使 得 
致 的 搜索 , 有 利于 


定 全 局 最 优 解 。 


= 0 x (X! — Xoo) (18) 


中 ，m 通常 取 值 为 0.01，X! 表 示 第 1 代 的 第 i 个 鸟巢 位 置 ， 


第 二 种 改进 方法 如 式 (19) 所 示 , 采用 了 文献 [18] 提 出 的 改 
方式 ， 主 要 通过 算法 的 迭代 次 数 对 步 长 因子 4 进行 控制 。 


在 
算 
搜 


算法 初期 ， 步 长 因子 4 的 值 较 大 ， 搜 索 范围 也 大 ， 降 低 了 
法 陷入 局 部 最 优 的 可 能 , 到 了 算法 后 期 , 步 长 自 适应 减 小 
索 精 度 得 到 了 提高 ， 更 利于 找到 最 优 解 : 


(19) 


fi 
C=(Qna 十 0 )xcos| 一 一 
L 


max 


™ 


为 [-0.05,0.05] 上 随机 步 长 因子 。 
2.4.2 信息 交流 

于 三 个 子 群 是 分 别 独立 的 进行 寻 优 操作 ， 为 了 提高 寻 
优 效率 ， 每 进化 上 代 后 ， 各 个 种 群 之 间 通 过 一 定 的 策略 进行 
种 群 间 个 体 的 信息 交流 ， 以 交换 种 群 信息 ,保持 种 群 的 多 样 性 。 
言 息 交 流 的 原则 是 : 将 子 群 中 的 较 差 个 体 与 最 优 个 体 进行 信 
息 交 流 , 从 而 引导 较 差 个 体 向 全 局 最 优 个 体位 置 进行 搜索 。 本 


其 中 ,4 表示 当前 迭代 次 数 , me 表示 连 代 总 次 数 , 取 ow =09 ， 


Step4 


优 鸟巢 。 按 照 式 (17) 进 行 
成 相同 数量 的 新 解 。 
每 兴 代 kk 次， 引入 差分 算 子 


的 信 居 交 Ly 用 改进 解码 算法 对 
更 新 当前 最 优 鸟 巢 。 


随机 将 种 群 分 为 3 个 


标准 Lévy 飞行 公式 和 两 种 改进 Lévy 飞 


鸟巢 位 置 ， 解 码 得 到 完工 时 间 ， 更 新 当前 最 


好 随机 游 走 , 淘汰 部 分 差 的 解 并 生 


来 实现 三 个 子 群 之 间 
其 进行 解码 得 到 完工 时 间 ， 


Step5 当 达 到 最 大 迭代 次 数 则 输出 全 局 最 优 解 ， 否 则 转 
到 Step3 继续 进行 迭代 。 


3 ”实例 仿真 与 分 析 


3.1 


实验 环境 和 测试 实例 


本 文采 用 Matlab R2018a 编程 ， 运 行 环 境 为 Intel(R) 
Xeon(R) CPU E5- @3.30GHz，RAM 8GB。 考 虑 到 目前 针对 


DCRFJSP 的 研究 成 
为 了 更 好 的 通过 比较 


果 较 少 ， 且 


实验 结果 来 验证 算法 的 性 和 


| 没有 标准 算 例 可 


Brandimarte[19 标 ; 准 测试 自 
测试 , 同时 根据 文献 [11] 中 的 了 
体 工 人 机 器 分 配 情况 如 


3.2 结果 和 分 析 


为 了 保证 算法 以 较 优 和 


的 状态 运行 


中 的 MK1-MK10 算 例 对 算法 进行 


供 参 考 比较 ， 
EE， 本 文选 择 


对 参数 进行 设 


的 正 交 实验 法 


文 算法 求解 MK01 算 
图 5(a) 为 最 大 迭代 次 数 iter 的 影 中 


例 ， 运 行 10 


[人 分 配方 式 进行 仿真 实验 , 具 
0 表 3 所 示 。 


由 于 算法 的 参数 对 算法 的 性 能 以 及 运行 时 间 有 很 大 影响 ， 
， 本 文 引入 文献 [20] 所 采用 
。 图 5 为 不 同 参数 情况 下 ， 本 
次 的 平均 值 变化 趋势 图 。 


二 


当 和 迭代 次 数 设 置 


较 小 时 ， 可 


台 E 扎 


能 导致 算法 在 尚未 收敛 的 情况 下 


被 迫 停止 , 从 而 影响 算法 的 寻 优 结果 , 随 着 迭代 次 数 的 增加 ， 


获得 的 解 的 质量 也 逐渐 提高 ,但 是 当 迭 代 次 数 到 


文 


引入 差分 进化 算法 的 DE/best/1 变异 策略 对 子 群 进行 差分 ， 


继续 增加 烛 代 次 数 


不 能 明显 提高 


增加 算法 的 复杂 度 ， 因 此 郁 
大 和 迭代 次 数 zer=200 。 图 5(b) 可 以 看 出 ， 当 种 群 


E 保 这 


线 图 , 从 图 中 可 以 看 出 ， 


了 茶 一 值 后 ， 


取得 最 优 解 的 次 数 ， 还 会 
运行 效率 的 基础 上 ， 设 置 最 


规模 在 50 左 
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右 ， 能 获得 比较 优秀 的 解 ， 因 此 设置 种 群 数量 =50 。 图 5(c) 
为 发 现 概率 Pe 的 影响 曲线 图 ， 当 发 现 概率 在 0.25 左右 时 ， 
得 到 的 解 最 优 ， 随 着 发 现 概 率 逐 渐 增 大 ， 得 到 的 解 的 质量 
逐渐 降低 ， 因 此 最 终 设置 发 现 概率 re=0.25 。 
表 3 工人 和 机 器 分 配 
Tab. 3 Worker and Machine Allocation 


改进 布谷 乌 算 法 求解 双 资 源 约束 柔性 车 间 调 度 问题 


Instance w Allocation 


M(W)={Mi, M3, Ms} , MW)={M,, Ma, Ms}, 
MW)={Mi, M4, Me} , M (Ws)={M,, Ms, Ma} 
MW)={Mi,Ms}, MW)={M,, M4}, 
MW)={M, Ma, Me} , MW)={M,, Ms, Me, M1}, 
M(Ws)={Me, Mi, Ms}, M(We)={M;s, Ms} 
MW)={Mi,M;3, M4}, MW)={M,, M4}, 
MW)={M,M,,M;} 
MW)={Mi, Ms, Mio} , MW)={M,, M1, M11}, 
MW)={M3, Ma, Mo, M1}, MWa)={M,, Mo, Mi, ， 
MWs)={Ms6, Mi, Ms, Mis}, MWe)={M;, Ms,Mio}, 
MW)={M, Mo, Ma, Mis}, MWe)= {Mi,M;, Mio, Mia} 
MW)={Mi,M;3, Ms}, M(W,)={M,,M,}, 
MW)={M;, M4}, MW)= {M1,M,, Ms} 
MW)={Mi,M;,Ms}, M(W)={M,, Ma,M,}, 
MW)={M3, M4, Ms, Mio} , MW) = {M1, M1,M,}, 
M(Ws)={Ms,Me, M37}, M(We)={M,, Ma, Ms, Mio} 

种 群 规模 n 


MK1-2 4 
MK3-4 6 


MK5 3 


MK6,10 8 


MK7,11 4 


MK 8,9,12, 
6 


夫 代 次 数 iter 


发 现 概率 pa 


02 03 04 05 06 07 08 09 1 
(0) 


图 5 参数 影响 趋势 图 
Fig.$5 Parameter influence trend graph 

为 测试 本 文 算法 的 有 效 性 ,将 其 与 不 同 的 算法 进行 了 
比较 ， 每 个 算 例 连续 充分 实验 20 次 ， 记 录 多 次 实验 的 结 
果 并 进行 分 析 , 仿真 实验 结果 见 表 4。 其 中 ，VNS 、KGFOA 
和 MBSA 分 别 为 文献 [10]、 文 献 [11] 和 文献 [21] 中 提出 的 算 
法 求解 DRCFJSP 所 得 的 最 优 值 ， 71Cs 为 普通 解码 方式 的 
改进 布谷 鸟 算法 求解 所 得 的 最 优 值 ， CSND 为 具有 改进 角 
码 方 式 的 标准 布谷 鸟 算法 求解 所 得 的 最 优 值 ，7CSND 为 
文采 用 的 具有 改进 解码 方式 的 改进 布谷 鸟 算法 求解 所 
的 最 优 值 。 

从 表 4 可 以 看 出 ,在 仅 对 布谷 鸟 算法 进行 改进 的 情况 下 
当 问 题 规模 较 小 时 ， 容 易 得 到 比较 优秀 的 值 ， 但 是 当 问 题 规 
模 逐 渐 增 大 ， 机 器 和 工人 空闲 时 间 也 随 之 增多 ， 普 通 的 解码 
方式 无 法 有 效 利 用 这 些 空 闲 时 间 ， 因 此 仅 靠 改进 布谷 岛 算 法 
已 经 无 法 得 到 更 好 的 解 。 对 于 CSND ， 由 于 布谷 岛 算 法 本 身 具 
有 一 定 的 优越 性 ， 能 较 好 地 兼顾 全 局 和 局 部 搜索 ， 再 加 上 改 
进 解码 方式 有 效 的 缩短 了 完工 时 间 ， 求 解 的 效率 也 得 到 了 一 
定 的 提升 。 对 于 本 文 所 使 用 的 ICSND ， 不 仅 对 布谷 鸟 算法 进 
行 了 改进 ， 增 强 了 算法 的 搜索 能 力 ， 还 对 解码 方式 进行 了 调 
整 ， 有 效 利 用 了 工人 和 机 器 的 空闲 时 间 ， 极 大 的 缩短 了 完工 
时 间 ， 将 ICSND 与 其 他 算法 进行 对 比 ， 可 以 发 现在 大 部 分 算 
例 的 求解 上 ， 7CSND 都 能 得 到 更 优 的 解 。 
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表 4 不 同 算 法 求解 Brandimarte 算 例 的 结果 对 比 


Tab.4 Comparison ofthe results of different algorithms for solving 


the Brandimarte example 


Instance VNS KGFOA MBSA4A ICS CSND JICSND 
MKI1 63 6061 47 49 46 44 
MK2 3 52 39 40 40 38 
MK3 204 204 204 213 204 204 
MK4 69 67 79 80 75 67 
MK5 337 287 239 244 237 235 
MK6 89 86 9 82 76 73 
MK7 184 184 185 198 185 181 
MK8 536 551 561 612 372 537 
MK9 437 407 407 457 436 430 
MK10 328 315 311 332 292 290 


图 6 和 图 7 为 采用 本 文 ICSND 算法 求解 算 例 MK1 和 MK2 
所 得 的 甘 特 图 ， 其 中 ， 横 坐标 为 加 工时 长 ， 纵 坐标 为 机 器 ， 
方 框 里 的 第 一 个 数字 表示 工件 ,第 二 个 数字 表示 对 应 的 工人 。 


MKD1 其 例 甘 竺 图 
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图 6 MKI 甘 特 图 
Fig.6 MK1 gantt chart 


MK02 因 例 甘 特 图 
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图 7 MK2 甘 特 图 
Fig.7 MK2 gantt chart 
为 了 更 好 的 衡量 本 文 所 使 用 的 算法 性 能 ， 引 入 文献 [10] 
所 使 用 的 相对 百分比 偏差 最 优 值 MRPP 和 相对 百分比 偏差 平 
均值 4RPD 这 两 个 指标 来 评估 算法 能 力 ， 计 算 公式 为 


miow 

MRPD=min C=C x100] (21) 
1 Cnrw ~ Ces 

ARPD = > Be 1o0] (22) 


其 中 ，Cw 为 各 个 算法 求 出 的 最 优 解 ，C% 为 本 文 求 出 的 最 优 解 。 

从 表 5 可 以 看 出 ， 采 用 了 改进 解码 方式 的 算法 MRPD 和 
ARPD 都 优 于 普通 解码 方式 的 算法 , 而 ICSND 无 论 从 解 的 质量 
还 是 最 优 解 的 获取 都 优 于 CSND , 改进 后 的 布谷 鸟 算 法 性 能 整 
体 优 于 标准 布谷 鸟 算 法 ， 更 强 的 搜索 能 力 。 表 6 通过 将 
ICSND 与 VVS 、KGFOA 和 MBSA 的 实验 结果 进行 对 比 可 以 看 出 ， 
ICSND 的 最 小 百分比 偏差 都 为 0， 并 且 平 均 百 分 比 偏差 也 更 
小 ， 得 到 的 解 更 加 稳定 ， 质 量 也 比较 好 。 而 其 他 三 种 算法 ， 
在 大 部 分 情况 下 的 求解 质量 都 不 如 ICSND 算法 ， 平 均 百分比 
偏差 的 值 也 更 大 。 因 此 相对 于 其 他 三 种 算法 ，7CSND 算法 具 


录用 定稿 


有 更 好 的 求解 能 力 ， 在 求解 DRCFJSP 问题 时 的 稳定 性 和 解 
的 质量 都 更 优 。 
表 5 不 同 布谷 鸟 算法 之 间 的 对 比分 析 


Tab. 53 Comparative analysis between different cuckoo algorithms 


ICSND 7TCSND 7TCS 


Instance Clow 
MRPD MRPD MRPD ARPD MRPD ARPD 


MKI1 44 0 :| 0 6.44 8.89 4.67 
MK2 38 0 pA 2.63 4.2 5.26 0.26 
MK3 204 0 2.45 0 3.34 5.39 0.69 
MK4 67 0 4.10 10.45 18.20 20.89 24.93 
MK5 235 0 8.89 1.70 2.46 4.68 6.09 
MK6 73 0 273 2.73 4.52 13.69 6.98 
MK7 181 0 2.70 2.21 3.20 6.07 9.66 
MK8 337 0 2.52 2.49 4.12 8.71 和 .72 
MK9 430 0 1.03 0.69 2.14 10.80 2.48 
MK10 290 0 0.89 0.34 1.44 .3 让 炒 3 


表 6 布谷 鸟 算法 和 其 他 三 种 算法 的 对 比分 析 
Tab.6 Comparative analysis of the cuckoo algorithm and 


the other three algorithms 


ICSND VNS KGFOA MBSA 


Instance MXNnxw Clow 


MRPD ARPD MRPD ARPD MRPD ARPD MRPD ARPD 


MKI1 10x6x4 44 0 3.11 35.55 58.76 40.00 46.67 0 5.11 
MK2 10x6x4 38 2.37 36.84 68.38 34.21 47.36 0 2.89 
MK3 1Sx8x6 204 2.45 0 49.18 0 41.00 0 4.49 
MK4 1Sx8x6 67 4.10 0 31.76 -11.94 26.28 5.94 15.52 
MK5 1Sx4x3 235 8.89 22.12 38.70 43.40 33.37 1.27 2.04 


MK6 10xlSx8 73 2.73 17.80 70.44 21.91 61.27 1.36 3.56 


© © © © © © 


MK7 20xSx4 181 .70 L665 9 5 _ 213 2 588 
MK8 20xl0x6 557 0 2.52 -1.95 20.79 -4.62 16.70 1.54 4.77 
MK9 20x10x6 430 0 1.03 -6.43 26.98 0.45 19.91 3.41 9.84 


MK10 20xlSx8 290 0 0.89 8.62 56.55 13.10 48.91 3.80 8.47 


4 ”结束 语 


本 文 在 布谷 鸟 算法 核心 框架 不 变 的 基础 下 对 其 进行 改进 
通过 将 布谷 鸟 种 群 划 分 为 3 个 子 群 ,并 分 别 采 用 不 同 Lévy 飞 
行 方式 进行 寻 优 , 同 时 引入 差分 算 子 进行 子 群 间 的 信息 交流 ， 
大 大 增加 了 算法 的 搜索 能 力 。 在 解码 时 设计 了 一 种 改进 解码 
方式 ， 在 避免 机 器 和 工人 的 加 工时 间 冲 突 的 同时 ， 主 动 寻找 
可 进行 插入 的 空闲 时 间 , 大 大 的 缩短 了 整个 加 工 的 完工 时 间 。 
最 后 通过 基准 测试 算 例 进行 实验 验证 和 比较 ， 证 明了 改进 布 
谷 鸟 算 法 和 改进 解码 方式 的 有 效 性 和 稳定 性 。 在 后 续 研 究 中 ， 
考虑 将 生产 过 程 中 常见 的 动态 因素 考虑 进去 , 比如 机 器 故障 、 
机 器 维护 、 紧 急 插 单 以 及 紧急 撤 单 等 。 
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